ASTIA  FILE  COPY 


a NEW  TEST  OF  COMPOUND  SYMMETRY 


by 


S.  N.  ROY 

INSTITUTE  OF  STATISTICS,  UNIVERSITY  OF  NORTH  CAROLINA 


Institute  of  Statistics 
Miseograph  Series  No.  97 

March  11,  1951* 


Work  sponsored  by  the  Office  of  Naval  Research  under  Contract  No. 
N7-onr-28h02. 


THIS  REPORT  HAS  BEEN  DELIMITED 
AND  CLEARED  FOR  PUBLIC  RELEASE 
UNDER  DOD  DIRECTIVE  5200.20  AND 
NO  RESTRICTIONS  ARE  IMPOSED  UPON 
ITS  USE  AND  DISCLOSURE. 

DISTRIBUTION  STATEMENT  A 

APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED,  . 


A NEW  TEST  OF  COMPOUND  SYMMETRY 


by 
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1.  Summary.  If  x^  and  x g have  a bivariate  normal  distribution  with 

a correlation  coefficient  p and  the  same  standard  deviations  a for  both, 
then  it  is  well  known  and  easy  to  check  that  x^  + Xg  and  x1  - Xg  are  un- 


correlated, which  forms  the  basis  of  Pitman’s  well-known  test  of  H^:  o = o^, 

for  a bivariate  normal  population,  in  terms  of  the  correlation  coefficient 
r between  x^  + Xg  a..d  x^  - x^  in  a random  sample  of  size,  say,  n,  from  this 


population.  Starting  from  this  test  which  has  a number  of  reasonably  good 
properties,  and  then  using  the  union- intersect ion  principle  £~1,  2 J , a 


test  is  obtained  for  compound  symmetry,  i.e.,  for  a ^ * o 


22 


PP 


and  all  o^'s  are  equal  (i  / j = 1,  2,  ...,  p),  where  is  any  element 


of  the  covariance  matrix  L of  a p-variate  normal  population. 

2.  Test  Construction.  Let  p(x^,  Xg)  denote  the  population  correlation 


coefficient  between  x,  and  Xg  and  r(x^,  Xg)  the  same  in  a random  sample 


of  size,  say  n,  from  that  population.  Then,  for  a set  of  stochastic 
variables  x'(l  x p)(=  x^,  ...,  x^)  having  a p-variate  normal  distribution, 

notice  that  for  any  arbitrary  non-null  vector  a’(l  x p),  a'(l  x p)  x(p  x 1) 

P P 

and  L x have  a bivariate  normal  distribution  and,  now  letting  L a.  = 0, 
i=l  1 i=l  1 
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P 

consider  the  hypothesis:  p(  L x , a'x)  = 0 = H (say).  Next  notice  that 

i=l  1 “ " - 


(2.1)  Hn:  ail  a 's  are  equal  and  all  o.  ^*s  are  equal  (i  -f  j = 1,  p) 


= AaHOa  = ^a  E P ( J^i'  * °J  > 

P 

where  a'(l  x p)  is  any  arbitrary  row  vector  subject  to  L a = 0. 

i=l  1 

Now  going  back  to  H , we  have  for  thiB  hypothesis  the  Pitman 
exotical  region,  say  W (a),  of  size  a,  given  by 

SL 

o p p 

(2.2)  W (a):  r“ ( L x , a'  x)  > r (n-2), 

- i=l  1 u 


where  r, (n-2)  is  the  upper  a/2-point  of  the  central  r -distribution  in 
random  samples  of  size  n. 

Hence,  by  the  union-intersection  heuristic  principle  £~1,  2_7 
we  have,  for  H ( = f\  H ),  the  critical  region  W(p)  of  Bize  p given 

u cl  US 

by 

(2.3)  V<(3):^~\a  /~r2(  L x , a'x)  > r^(n-2)_7  , 

- i=l 

2 P p 

i.e,,  Supa  r ( h x^,  a'x)  > r,  (n-2) 


P P P-1  P 

i.e.,  Sup  r ( Lx,  L a (x  - x ) ) > r^(n-2)  , 
- i=l  1 i=l  1 1 P C 

P P-1 

remembering  that  L a.  =0,  i.e., a =-L  a. 

i=l  1 P i=l  1 
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It  Is  easy  to  check  obat 


(2.4) 


P-1 


Sup  r ( L x , £ a (x  - x ) ) 

- i=l  1 1=1  1 1 p 


= square  of  the  sample  multiple  correlation  between  L x,  and  the  (p  - i)- 

i=l  x 

Set  Of  vnrjahles  (x.  - X ).  (x.  - X ),  ....  (x  ,-x)= 

' 1 p’  ' • * p ' p-1  p' 

2 P 2 

R C L \ and  (Xx  - xp)j  (x2  - xp),  . ..,  (xpl  - xp)  J '■  R (say). 

Notice  that,  since  the  new  p-set  also  has  the  p-variate  noraal  distribution, 
this  R has  the  well-known  multiple  correlation  distribution  wvLth  degrees  of 
freedom  p-1  and  n - p and  a non-centrality  parameter  whies  . Lb  the  popula- 

P 

tion  multiple  correlation  between  1 x and  the  (p  - 1) -set;  above  and 

i=l 

2 

which  let  ue  call  p . It  is  easy  to  check  that  p = 0,  i.e,, , that  R has 
the  central  multiple  correlation  distribution,  if  and  onlyilf 
P 

p(  I x , x - x ) = 0 (i  = 1,  2,  ...,  p - 1),  i.e.,  if  ani  only  if 
i=l  1 1 p 


F P s 

p ( 1 x , a'x)  = 0 (for  all  non-null  a subject  to  h a --  0). 
i=l  ” i=-l  1 


We  have  thus,  for  testing  compound  symmetry,  the  crivtical  region 
W(P)  of  size  P given  by 
P 

(2.5)  R L L and  (x,  - xp),  ...,  (xp_1  - ,cp)  J > (p  - 1,  n - p), 
when  Rp  is  the  upper  P-point  **  the  well-known  central  mulOtiple  correlation 


distribution . 
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It  can  be  checked  after  some  little  algebra  that  in  terms,  respectively, 
of  the  elements  of  the  sample  and  population  covariance  matrices  S and  £, 

R and  p will  be  given  by 


(2.6) 


when 


R = 1 - o L 


, P P 1 

z.z1/  ( R z.)  ( L z ), 


i=l  i=l  i=l 

P ^ . P P 


P2 *  = 1 - P ^ S U ( £ C1)  , 

i=l  i=l  i=l 


p i p ii 

zi  * A 8i y 2 * A * - 


J-l 


J=1 


P i P 

t = L a and  £ = h 0 

* j i J x ^ 


ij 


J=1 


J=1 


Note  that  s . , * s . . , 
ij  Ji  * 


M = J1 


iJ  _ 


sJ  , o.  . * o,,  and  o **  * o' 
' 1J  Ji 


The  power  properties  of  this  test  will  be  discussed  in  a later  note. 
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